Blood DNA methylation profiling identifies cathepsin Z dysregulation in pulmonary arterial hypertension.

Pulmonary arterial hypertension (PAH) is characterised by pulmonary vascular remodelling causing premature death from right heart failure. Established DNA variants influence PAH risk, but susceptibility from epigenetic changes is unknown. We addressed this through epigenome-wide association study (EWAS), testing 865,848 CpG sites for association with PAH in 429 individuals with PAH and 1226 controls. Three loci, at Cathepsin Z (CTSZ, cg04917472), Conserved oligomeric Golgi complex 6 (COG6, cg27396197), and Zinc Finger Protein 678 (ZNF678, cg03144189), reached epigenome-wide significance (p < 10-7) and are hypermethylated in PAH, including in individuals with PAH at 1-year follow-up. Of 16 established PAH genes, only cg10976975 in BMP10 shows hypermethylation in PAH. Hypermethylation at CTSZ is associated with decreased blood cathepsin Z mRNA levels. Knockdown of CTSZ expression in human pulmonary artery endothelial cells increases caspase-3/7 activity (p < 10-4). DNA methylation profiles are altered in PAH, exemplified by the pulmonary endothelial function modifier CTSZ, encoding protease cathepsin Z.

Rare variant association analysis in 51,256 type 2 diabetes cases and 370,487 controls informs the spectrum of pathogenicity of monogenic diabetes genes.

We meta-analyzed array data imputed with the TOPMed reference panel and whole-genome sequence (WGS) datasets and performed the largest, rare variant (minor allele frequency as low as 5×10-5) GWAS meta-analysis of type 2 diabetes (T2D) comprising 51,256 cases and 370,487 controls. We identified 52 novel variants at genome-wide significance (p<5 × 10-8), including 8 novel variants that were either rare or ancestry-specific. Among them, we identified a rare missense variant in HNF4A p.Arg114Trp (OR=8.2, 95% confidence interval [CI]=4.6-14.0, p = 1.08×10-13), previously reported as a variant implicated in Maturity Onset Diabetes of the Young (MODY) with incomplete penetrance. We demonstrated that the diabetes risk in carriers of this variant was modulated by a T2D common variant polygenic risk score (cvPRS) (carriers in the top PRS tertile [OR=18.3, 95%CI=7.2-46.9, p=1.2×10-9] vs carriers in the bottom PRS tertile [OR=2.6, 95% CI=0.97-7.09, p = 0.06]. Association results identified eight variants of intermediate penetrance (OR>5) in monogenic diabetes (MD), which in aggregate as a rare variant PRS were associated with T2D in an independent WGS dataset (OR=4.7, 95% CI=1.86-11.77], p = 0.001). Our data also provided support evidence for 21% of the variants reported in ClinVar in these MD genes as benign based on lack of association with T2D. Our work provides a framework for using rare variant imputation and WGS analyses in large-scale population-based association studies to identify large-effect rare variants and provide evidence for informing variant pathogenicity.

Pre-RNA splicing in metabolic homeostasis and liver disease.

The liver plays a key role in sensing nutritional and hormonal inputs to maintain metabolic homeostasis. Recent studies into pre-mRNA splicing and alternative splicing (AS) and their effects on gene expression have revealed considerable transcriptional complexity in the liver, both in health and disease. While the contribution of these mechanisms to cell and tissue identity is widely accepted, their role in physiological and pathological contexts within tissues is just beginning to be appreciated. In this review, we showcase recent studies on the splicing and AS of key genes in metabolic pathways in the liver, the effect of metabolic signals on the spliceosome, and therapeutic intervention points based on RNA splicing.

RISING STARS: Liver sinusoidal endothelial transcription factors in metabolic homeostasis and disease.

Liver sinusoidal endothelial cells (LSECs) are highly specialised endothelial cells that form the liver microvasculature. LSECs maintain liver homeostasis, scavenging bloodborne molecules, regulating immune response, and actively promoting hepatic stellate cell quiescence. These diverse functions are underpinned by a suite of unique phenotypical attributes distinct from other blood vessels. In recent years, studies have begun to reveal the specific contributions of LSECs to liver metabolic homeostasis and how LSEC dysfunction associates with disease aetiology. This has been particularly evident in the context of non-alcoholic fatty liver disease (NAFLD), the hepatic manifestation of metabolic syndrome, which is associated with the loss of key LSEC phenotypical characteristics and molecular identity. Comparative transcriptome studies of LSECs and other endothelial cells, together with rodent knockout models, have revealed that loss of LSEC identity through disruption of core transcription factor activity leads to impaired metabolic homeostasis and to hallmarks of liver disease. This review explores the current knowledge of LSEC transcription factors, covering their roles in LSEC development and maintenance of key phenotypic features, which, when disturbed, lead to loss of liver metabolic homeostasis and promote features of chronic liver diseases, such as non-alcoholic liver disease.

Integrative genomic analyses in adipocytes implicate DNA methylation in human obesity and diabetes.

DNA methylation variations are prevalent in human obesity but evidence of a causative role in disease pathogenesis is limited. Here, we combine epigenome-wide association and integrative genomics to investigate the impact of adipocyte DNA methylation variations in human obesity. We discover extensive DNA methylation changes that are robustly associated with obesity (N = 190 samples, 691 loci in subcutaneous and 173 loci in visceral adipocytes, P < 1 × 10-7). We connect obesity-associated methylation variations to transcriptomic changes at >500 target genes, and identify putative methylation-transcription factor interactions. Through Mendelian Randomisation, we infer causal effects of methylation on obesity and obesity-induced metabolic disturbances at 59 independent loci. Targeted methylation sequencing, CRISPR-activation and gene silencing in adipocytes, further identifies regional methylation variations, underlying regulatory elements and novel cellular metabolic effects. Our results indicate DNA methylation is an important determinant of human obesity and its metabolic complications, and reveal mechanisms through which altered methylation may impact adipocyte functions.

Dax1 modulates ERα-dependent hypothalamic estrogen sensing in female mice.

Coupling the release of pituitary hormones to the developmental stage of the oocyte is essential for female fertility. It requires estrogen to restrain kisspeptin (KISS1)-neuron pulsatility in the arcuate hypothalamic nucleus, while also exerting a surge-like effect on KISS1-neuron activity in the AVPV hypothalamic nucleus. However, a mechanistic basis for this region-specific effect has remained elusive. Our genomic analysis in female mice demonstrate that some processes, such as restraint of KISS1-neuron activity in the arcuate nucleus, may be explained by region-specific estrogen receptor alpha (ERα) DNA binding at gene regulatory regions. Furthermore, we find that the Kiss1-locus is uniquely regulated in these hypothalamic nuclei, and that the nuclear receptor co-repressor NR0B1 (DAX1) restrains its transcription specifically in the arcuate nucleus. These studies provide mechanistic insight into how ERα may control the KISS1-neuron, and Kiss1 gene expression, to couple gonadotropin release to the developmental stage of the oocyte.

SOX17 Enhancer Variants Disrupt Transcription Factor Binding And Enhancer Inactivity Drives Pulmonary Hypertension.

BACKGROUND: Pulmonary arterial hypertension (PAH) is a rare disease characterized by remodeling of the pulmonary arteries, increased vascular resistance, and right-sided heart failure. Genome-wide association studies of idiopathic/heritable PAH established novel genetic risk variants, including conserved enhancers upstream of transcription factor (TF) SOX17 containing 2 independent signals. SOX17 is an important TF in embryonic development and in the homeostasis of pulmonary artery endothelial cells (hPAEC) in the adult. Rare pathogenic mutations in SOX17 cause heritable PAH. We hypothesized that PAH risk alleles in an enhancer region impair TF-binding upstream of SOX17, which in turn reduces SOX17 expression and contributes to disturbed endothelial cell function and PAH development. METHODS: CRISPR manipulation and siRNA were used to modulate SOX17 expression. Electromobility shift assays were used to confirm in silico-predicted TF differential binding to the SOX17 variants. Functional assays in hPAECs were used to establish the biological consequences of SOX17 loss. In silico analysis with the connectivity map was used to predict compounds that rescue disturbed SOX17 signaling. Mice with deletion of the SOX17-signal 1 enhancer region (SOX17-4593/enhKO) were phenotyped in response to chronic hypoxia and SU5416/hypoxia. RESULTS: CRISPR inhibition of SOX17-signal 2 and deletion of SOX17-signal 1 specifically decreased SOX17 expression. Electromobility shift assays demonstrated differential binding of hPAEC nuclear proteins to the risk and nonrisk alleles from both SOX17 signals. Candidate TFs HOXA5 and ROR-α were identified through in silico analysis and antibody electromobility shift assays. Analysis of the hPAEC transcriptomes revealed alteration of PAH-relevant pathways on SOX17 silencing, including extracellular matrix regulation. SOX17 silencing in hPAECs resulted in increased apoptosis, proliferation, and disturbance of barrier function. With the use of the connectivity map, compounds were identified that reversed the SOX17-dysfunction transcriptomic signatures in hPAECs. SOX17 enhancer knockout in mice reduced lung SOX17 expression, resulting in more severe pulmonary vascular leak and hypoxia or SU5416/hypoxia-induced pulmonary hypertension. CONCLUSIONS: Common PAH risk variants upstream of the SOX17 promoter reduce endothelial SOX17 expression, at least in part, through differential binding of HOXA5 and ROR-α. Reduced SOX17 expression results in disturbed hPAEC function and PAH. Existing drug compounds can reverse the disturbed SOX17 pulmonary endothelial transcriptomic signature.

Regulatory de novo mutations underlying intellectual disability.

The genetic aetiology of a major fraction of patients with intellectual disability (ID) remains unknown. De novo mutations (DNMs) in protein-coding genes explain up to 40% of cases, but the potential role of regulatory DNMs is still poorly understood. We sequenced 63 whole genomes from 21 ID probands and their unaffected parents. In addition, we analysed 30 previously sequenced genomes from exome-negative ID probands. We found that regulatory DNMs were selectively enriched in fetal brain-specific enhancers as compared with adult brain enhancers. DNM-containing enhancers were associated with genes that show preferential expression in the prefrontal cortex. Furthermore, we identified recurrently mutated enhancer clusters that regulate genes involved in nervous system development (CSMD1, OLFM1, and POU3F3). Most of the DNMs from ID probands showed allele-specific enhancer activity when tested using luciferase assay. Using CRISPR-mediated mutation and editing of epigenomic marks, we show that DNMs at regulatory elements affect the expression of putative target genes. Our results, therefore, provide new evidence to indicate that DNMs in fetal brain-specific enhancers play an essential role in the aetiology of ID.

Liver RBFOX2 regulates cholesterol homeostasis via Scarb1 alternative splicing in mice.

RNA alternative splicing (AS) expands the regulatory potential of eukaryotic genomes. The mechanisms regulating liver-specific AS profiles and their contribution to liver function are poorly understood. Here, we identify a key role for the splicing factor RNA-binding Fox protein 2 (RBFOX2) in maintaining cholesterol homeostasis in a lipogenic environment in the liver. Using enhanced individual-nucleotide-resolution ultra-violet cross-linking and immunoprecipitation, we identify physiologically relevant targets of RBFOX2 in mouse liver, including the scavenger receptor class B type I (Scarb1). RBFOX2 function is decreased in the liver in diet-induced obesity, causing a Scarb1 isoform switch and alteration of hepatocyte lipid homeostasis. Our findings demonstrate that specific AS programmes actively maintain liver physiology, and underlie the lipotoxic effects of obesogenic diets when dysregulated. Splice-switching oligonucleotides targeting this network alleviate obesity-induced inflammation in the liver and promote an anti-atherogenic lipoprotein profile in the blood, underscoring the potential of isoform-specific RNA therapeutics for treating metabolism-associated diseases.

GWAS meta-analysis of intrahepatic cholestasis of pregnancy implicates multiple hepatic genes and regulatory elements.

Intrahepatic cholestasis of pregnancy (ICP) is a pregnancy-specific liver disorder affecting 0.5-2% of pregnancies. The majority of cases present in the third trimester with pruritus, elevated serum bile acids and abnormal serum liver tests. ICP is associated with an increased risk of adverse outcomes, including spontaneous preterm birth and stillbirth. Whilst rare mutations affecting hepatobiliary transporters contribute to the aetiology of ICP, the role of common genetic variation in ICP has not been systematically characterised to date. Here, we perform genome-wide association studies (GWAS) and meta-analyses for ICP across three studies including 1138 cases and 153,642 controls. Eleven loci achieve genome-wide significance and have been further investigated and fine-mapped using functional genomics approaches. Our results pinpoint common sequence variation in liver-enriched genes and liver-specific cis-regulatory elements as contributing mechanisms to ICP susceptibility.

Dysregulated RNA polyadenylation contributes to metabolic impairment in non-alcoholic fatty liver disease.

Pre-mRNA processing is an essential mechanism for the generation of mature mRNA and the regulation of gene expression in eukaryotic cells. While defects in pre-mRNA processing have been implicated in a number of diseases their involvement in metabolic pathologies is still unclear. Here, we show that both alternative splicing and alternative polyadenylation, two major steps in pre-mRNA processing, are significantly altered in non-alcoholic fatty liver disease (NAFLD). Moreover, we find that Serine and Arginine Rich Splicing Factor 10 (SRSF10) binding is enriched adjacent to consensus polyadenylation motifs and its expression is significantly decreased in NAFLD, suggesting a role mediating pre-mRNA dysregulation in this condition. Consistently, inactivation of SRSF10 in mouse and human hepatocytes in vitro, and in mouse liver in vivo, was found to dysregulate polyadenylation of key metabolic genes such as peroxisome proliferator-activated receptor alpha (PPARA) and exacerbate diet-induced metabolic dysfunction. Collectively our work implicates dysregulated pre-mRNA polyadenylation in obesity-induced liver disease and uncovers a novel role for SRSF10 in this process.

Deletion of Regulatory Elements with All-in-One CRISPR-Cas9 Vectors.

Loss-of-function experiments are essential for the functional investigation of cis-regulatory elements (CREs), such as transcriptional enhancers. This can be achieved with CRISPR-Cas9 using pairs of single guide RNAs (sgRNAs) to target the flanking regions of a CRE. Here, I describe a single-step protocol to rapidly and inexpensively generate vectors co-expressing two sgRNAs, which allows re-usage of gRNAs oligonucleotides from one experimental design to another. This protocol is applicable to cloning sgRNAs into virtually any CRISPR-Cas9 backbone that allows cloning using Golden Gate, by adapting the primer design.

Epigenetics of Hepatic Insulin Resistance.

Insulin resistance (IR) is largely recognized as a unifying feature that underlies metabolic dysfunction. Both lifestyle and genetic factors contribute to IR. Work from recent years has demonstrated that the epigenome may constitute an interface where different signals may converge to promote IR gene expression programs. Here, we review the current knowledge of the role of epigenetics in hepatic IR, focusing on the roles of DNA methylation and histone post-translational modifications. We discuss the broad epigenetic changes observed in the insulin resistant liver and its associated pathophysiological states and leverage on the wealth of 'omics' studies performed to discuss efforts in pinpointing specific loci that are disrupted by these changes. We envision that future studies, with increased genomic resolution and larger cohorts, will further the identification of biomarkers of early onset hepatic IR and assist the development of targeted interventions. Furthermore, there is growing evidence to suggest that persistent epigenetic marks may be acquired over prolonged exposure to disease or deleterious exposures, highlighting the need for preventative medicine and long-term lifestyle adjustments to avoid irreversible or long-term alterations in gene expression.

Chromatin 3D interaction analysis of the STARD10 locus unveils FCHSD2 as a regulator of insulin secretion.

Using chromatin conformation capture, we show that an enhancer cluster in the STARD10 type 2 diabetes (T2D) locus forms a defined 3-dimensional (3D) chromatin domain. A 4.1-kb region within this locus, carrying 5 T2D-associated variants, physically interacts with CTCF-binding regions and with an enhancer possessing strong transcriptional activity. Analysis of human islet 3D chromatin interaction maps identifies the FCHSD2 gene as an additional target of the enhancer cluster. CRISPR-Cas9-mediated deletion of the variant region, or of the associated enhancer, from human pancreas-derived EndoC-ßH1 cells impairs glucose-stimulated insulin secretion. Expression of both STARD10 and FCHSD2 is reduced in cells harboring CRISPR deletions, and lower expression of STARD10 and FCHSD2 is associated, the latter nominally, with the possession of risk variant alleles in human islets. Finally, CRISPR-Cas9-mediated loss of STARD10 or FCHSD2, but not ARAP1, impairs regulated insulin secretion. Thus, multiple genes at the STARD10 locus influence ß cell function.

Liver gene regulatory networks: Contributing factors to nonalcoholic fatty liver disease.

Metabolic diseases such as nonalcoholic fatty liver disease (NAFLD) result from complex interactions between intrinsic and extrinsic factors, including genetics and exposure to obesogenic environments. These risk factors converge in aberrant gene expression patterns in the liver, which are underlined by altered cis-regulatory networks. In homeostasis and in disease states, liver cis-regulatory networks are established by coordinated action of liver-enriched transcription factors (TFs), which define enhancer landscapes, activating broad gene programs with spatiotemporal resolution. Recent advances in DNA sequencing have dramatically expanded our ability to map active transcripts, enhancers and TF cistromes, and to define the 3D chromatin topology that contains these elements. Deployment of these technologies has allowed investigation of the molecular processes that regulate liver development and metabolic homeostasis. Moreover, genomic studies of NAFLD patients and NAFLD models have demonstrated that the liver undergoes pervasive regulatory rewiring in NAFLD, which is reflected by aberrant gene expression profiles. We have therefore achieved an unprecedented level of detail in the understanding of liver cis-regulatory networks, particularly in physiological conditions. Future studies should aim to map active regulatory elements with added levels of resolution, addressing how the chromatin landscapes of different cell lineages contribute to and are altered in NAFLD and NAFLD-associated metabolic states. Such efforts would provide additional clues into the molecular factors that trigger this disease. This article is categorized under: Biological Mechanisms > Metabolism Biological Mechanisms > Regulatory Biology Laboratory Methods and Technologies > Genetic/Genomic Methods.

Pancreatic Islet Transcriptional Enhancers and Diabetes.

PURPOSE OF REVIEW: Common genetic variants that associate with type 2 diabetes risk are markedly enriched in pancreatic islet transcriptional enhancers. This review discusses current advances in the annotation of islet enhancer variants and their target genes. RECENT FINDINGS: Recent methodological advances now allow genetic and functional mapping of diabetes causal variants at unprecedented resolution. Mapping of enhancer-promoter interactions in human islets has provided a unique appreciation of the complexity of islet gene regulatory processes and enabled direct association of noncoding diabetes risk variants to their target genes. The recently improved human islet enhancer annotations constitute a framework for the interpretation of diabetes genetic signals in the context of pancreatic islet gene regulation. In the future, integration of existing and yet to come regulatory maps with genetic fine-mapping efforts and in-depth functional characterization will foster the discovery of novel diabetes molecular risk mechanisms.

Human pancreatic islet three-dimensional chromatin architecture provides insights into the genetics of type 2 diabetes.

Genetic studies promise to provide insight into the molecular mechanisms underlying type 2 diabetes (T2D). Variants associated with T2D are often located in tissue-specific enhancer clusters or super-enhancers. So far, such domains have been defined through clustering of enhancers in linear genome maps rather than in three-dimensional (3D) space. Furthermore, their target genes are often unknown. We have created promoter capture Hi-C maps in human pancreatic islets. This linked diabetes-associated enhancers to their target genes, often located hundreds of kilobases away. It also revealed >1,300 groups of islet enhancers, super-enhancers and active promoters that form 3D hubs, some of which show coordinated glucose-dependent activity. We demonstrate that genetic variation in hubs impacts insulin secretion heritability, and show that hub annotations can be used for polygenic scores that predict T2D risk driven by islet regulatory variants. Human islet 3D chromatin architecture, therefore, provides a framework for interpretation of T2D genome-wide association study (GWAS) signals.

Genetic determinants of risk in pulmonary arterial hypertension: international genome-wide association studies and meta-analysis.

BACKGROUND: Rare genetic variants cause pulmonary arterial hypertension, but the contribution of common genetic variation to disease risk and natural history is poorly characterised. We tested for genome-wide association for pulmonary arterial hypertension in large international cohorts and assessed the contribution of associated regions to outcomes. METHODS: We did two separate genome-wide association studies (GWAS) and a meta-analysis of pulmonary arterial hypertension. These GWAS used data from four international case-control studies across 11 744 individuals with European ancestry (including 2085 patients). One GWAS used genotypes from 5895 whole-genome sequences and the other GWAS used genotyping array data from an additional 5849 individuals. Cross-validation of loci reaching genome-wide significance was sought by meta-analysis. Conditional analysis corrected for the most significant variants at each locus was used to resolve signals for multiple associations. We functionally annotated associated variants and tested associations with duration of survival. All-cause mortality was the primary endpoint in survival analyses. FINDINGS: A locus near SOX17 (rs10103692, odds ratio 1·80 [95% CI 1·55-2·08], p=5·13 × 10-15) and a second locus in HLA-DPA1 and HLA-DPB1 (collectively referred to as HLA-DPA1/DPB1 here; rs2856830, 1·56 [1·42-1·71], p=7·65 × 10-20) within the class II MHC region were associated with pulmonary arterial hypertension. The SOX17 locus had two independent signals associated with pulmonary arterial hypertension (rs13266183, 1·36 [1·25-1·48], p=1·69 × 10-12; and rs10103692). Functional and epigenomic data indicate that the risk variants near SOX17 alter gene regulation via an enhancer active in endothelial cells. Pulmonary arterial hypertension risk variants determined haplotype-specific enhancer activity, and CRISPR-mediated inhibition of the enhancer reduced SOX17 expression. The HLA-DPA1/DPB1 rs2856830 genotype was strongly associated with survival. Median survival from diagnosis in patients with pulmonary arterial hypertension with the C/C homozygous genotype was double (13·50 years [95% CI 12·07 to >13·50]) that of those with the T/T genotype (6·97 years [6·02-8·05]), despite similar baseline disease severity. INTERPRETATION: This is the first study to report that common genetic variation at loci in an enhancer near SOX17 and in HLA-DPA1/DPB1 is associated with pulmonary arterial hypertension. Impairment of SOX17 function might be more common in pulmonary arterial hypertension than suggested by rare mutations in SOX17. Further studies are needed to confirm the association between HLA typing or rs2856830 genotyping and survival, and to determine whether HLA typing or rs2856830 genotyping improves risk stratification in clinical practice or trials. FUNDING: UK NIHR, BHF, UK MRC, Dinosaur Trust, NIH/NHLBI, ERS, EMBO, Wellcome Trust, EU, AHA, ACClinPharm, Netherlands CVRI, Dutch Heart Foundation, Dutch Federation of UMC, Netherlands OHRD and RNAS, German DFG, German BMBF, APH Paris, INSERM, Université Paris-Sud, and French ANR.